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\ Abstract. This is the first in a series of papers on scattering theory for one- 

• dimensional Schrodinger operators with highly singular potentials q G H{^l(R). In 

, this paper, we study Miura potentials q associated to positive Schrodinger operators 

that admit a Riccati representation q = u' + for a unique u e L^(R) n L^(R). 
Such potentials have a well-defined reflection coefhcient r{k) that satisfies |r(fc)| < 1 
and determines u uniquely. We show that the scattering map S : u 1— > r is real- 
analytic with real-analytic inverse. To do so, we exploit a natural complexification 
^ . of the scattering map associated with the ZS-AKNS system. In subsequent papers, 

^ I we will consider larger classes of potentials including singular potentials with bound 

states. 
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1. Introduction 

This is the first of a series of papers on scattering theory for one-dimensional Schrodinger 
operators 

S:= + ^ 
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on the line with highly singular potentials q. Our goal is to extend the inverse scattering 
method as discussed e.g. in [13,17,29-31] in order to study initial value problems for 
completely integrable dispersive equations with highly singular initial data. 

It is natural to begin with the case where the Schrodinger operator has no bound 
states; the corresponding potentials then admit a Riccati representation given by the 
Miura map [32]. Recall that the Miura map is the nonlinear mapping 

u ^ u + u . 

It is not difficult to see [32] that if u{x,t) is a smooth solution of the mKdV equation, 
then 

du 

(Bu)(x, t) = T— (x, t) + u^ix, t) 
ox 

is a smooth solution of the KdV equation. For this reason, the Miura map has played 
a fundamental role in the study of existence and well-posedness questions for these two 
equations. 

The range of the Miura map may be characterized as follows. For q G H^l(K) 
real- valued and G C^(R), consider the Schrodinger form 

5^:= I \v'{x)\'dx+{q,\^\'), (1.2) 

where ( ■ , ■ ) is the pairing between H^^l{R) and H^^^p{M). In [26] it was shown that 
if q is any real- valued distribution in H^l(R) for which the Schrodinger form s is 
nonnegative, then q may be presented as g = Bu for a function u G L^q^(]R) that 
need not be unique. We will call such a potential q a Miura potential, and we will call 
any function u with q = Bu a Riccati representative for q. It is not difficult to see that 
two Riccati representatives for a given q differ by a continuous function. Any Riccati 
representative u is the logarithmic derivative of a positive distributional solution to the 
zero-energy Schrodinger equation for q. 

In this paper, we will study inverse scattering for the subset of the Miura potentials 
which admit a Riccati representative u G L^(M) n L^(M). In this case, such a 
Riccati representative is unique, so we may parameterize the potentials by their Riccati 
representatives. In Paper II of this series [24], we consider Miura potentials q which 
do not have a unique Miura representative but instead possess the following property: 
There exist Riccati representatives m+ and which belong to L^(]R) and, in addition, 
is integrable on (— oo,0), while is integrable on (0,oo). It will turn out that a 
potential q in this class is uniquely determined by the data M-|(_ooo)' ^+I(ooo)' 
"jump" (m+ — u_)(0). This class includes the Faddeev-Marchenko class of potentials 
q G L^(]R, (1 + |a;|)dx) as well as many highly oscillatory and distributional potentials. 
Finally, in Paper III [25] , we will show how to add bound states. The class of potentials 
covered in the present paper is "non-generic" in the sense that, even for regular short- 
range potentials, the Riccati representative is generically non-unique; on the other hand. 



I See Subsection 12.11 for explicit definitions of all the function spaces appearing in the paper 
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the ideas used here lay the groundwork for the analysis of "generic" singular potentials 
in Papers II and III. 

The inverse scattering problem for Miura potentials with u G L^(]R) fl L^(]R) was 
worked out in detail from the point of view of Schrodinger scattering in the thesis of 
the first author [20]. In this paper, we take a somewhat different point of view and 
exploit a natural complexification of the problem that leads to the ZS-AKNS system 
(see Zakharov and Shabat [41] and Ablowitz, Kaup, Newell, and Segur [1]) associated to 
the defocusing nonlinear Schrodinger equation. Elsewhere [10], we will apply the refined 
continuity results obtained here to study solutions to the defocusing NLS equation with 
singular initial data. 

Let us first describe the connection between the ZS-AKNS system and scattering 
with Miura potentials. For a real- valued u in the Schwartz class S'(M), suppose that \l/ 
is a matrix-valued solution of the ZS-AKNS systen|§| 

m' = iza^ + q{x)^ , (1.3) 

where 

2 ^ \^ -\ ) \ u{x) 

An easy calculation shows that if \E' solves (11. 3p for a real- valued u G 5'(]R), then the 
row vector 

X:=(l 1)^ 
solves the Schrodinger equation 

-x" + <i{x)x = -jx, 

where 

q{x) = u'{x) + u^{x) (1.4) 

is the image of u under the Miura map (11. II) . Hence, in order to study the direct and 
inverse scattering maps for the Schrodinger equation with a Miura potential, it suffices 
to study the analogous maps for the ZS-AKNS system. 

To state the connection between the scattering maps more precisely, suppose again 
that u G S'(M) is real- valued and s G M. As is well-known (see for example [5], where 
more general problems under weaker regularity and decay assumptions are studied), 
there exist solutions \E'± of the ZS-AKNS system with respective asymptotics 

lim |\E'±(a;, s) — exp(isa;cr)| = (1.5) 

x— ►itoo 

§ Our equation (|1.3p is actually a special case of the ZS-AKNS system since in general one takes 

V ""2(2;) 

for complex- valued ui and U2- 
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and a unique matrix A(s) of the form 
A(.) = 



a{s) h{s) 
b{s) a{s) 



with 



such that 



\ais)\' - mi' 



^+(a:,s) = ^_(x,s)A(s). (1.6) 

If -u G 5'(]R), we have a — 1 G S{M.) and b G S'(M). The function a{s) extends to a 
bounded analytic function a{z) on Imz > with \a{z)\ ^ 1 as \z\ oo. 
The left reflection coefficient is given by 

r_{s) :=b{s)/a{s), (1.7) 

and |r_(s)| < 1 strictly if m G ^(IR). Similarly, the right reflection coefficient is given by 

r+(s) := -b{s)/a{s). (1.8) 

These reflection coefficients are identical (up to the change s ^ 2s in the spectral 
parameter) to the reflection coefficients in the Schrodinger problem for the Miura 
potential (11.41) . Indeed, if 

X± := (1 1)^± 

and x± = {xf^xf)^ we can use (11.51) and (11.61) to conclude that 

a(s)e'^^/2 ^ 5(s)e-iW2^ ^ -oo, 

X +00, 

and 



Xn^,s)r.{ -y^'^ ' ' ^^7-' (1.9) 



^^'^'^'"^ ^ ^ a(s)e--/2 _5(^)e-/2^ ^ _ +00.' ^^'^^^ 

Thus, up to the reparametrization, r_.(s) is the usual Schrodinger reflection coefficient 
for scattering from the left, and r+(s) is the usual Schrodinger reflection coefficient for 
scattering from the right. The asymptotic relations (11.91) and (11.101) show that, up to 
reparametrization, Xi{'^s) and X2i.'i^) the usual Jost solutions at +00 and —00 
respectively. 

These calculations motivate us to consider the inverse problem for the ZS-AKNS 
system, now for u G L^(R) fl L^(R). This class of u is chosen so that the corresponding 
Miura potentials lie in if~^(M) and define self-adjoint Schrodinger operators S [23,26]; 
moreover, they have sufficient decay to allow a well-defined scattering theory. It will 
pose no essential difficulty in considering the inverse scattering problem for the system 
fOD with 



Q(-) = \ r (1-11) 
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and u complex-valued; the scattering map for the Schrodinger problem with Miura 
potential q = u' + and u real- valued will be the restriction of the resulting scattering 
maps to real- valued u. 

For this more general scattering problem, there still exist matrix-valued Jost 
solutions and a matrix A with the above properties. In particular, the reflection 
coefficients r± are well defined, and we introduce the direct scattering maps S± via 

: u ^ r_. 

Define the Fourier transforms ^+ and !F-. via 



J — oo 



' — oo 

it will be convenient to work with 



:= (^+r+)(x) = ^ I e^^V(s) ds, (1.12) 



271 



-oo 
oo 



:= (J^_r_)(x) = — / e-'"^r_(s)ds (1.13) 

and to analyze the mappings u ^ and u ^ given by o S+ and JF_ o S- . 
To state our results, let us define the Banach space 

X := L^(M) nL2(M) (1.14) 
with the norm 

\\F\\x:=\\F\\l^ + \\F\\l2. (1.15) 

Let 

Xi := |fgX : ||F||oo < l} 
with norm induced from X, where 

/oo 
e-^^^F(a;) da; = {r~^F){s). (1.16) 
'OO 

Our first result is: 

Theorem 1.1 The mappings o 5+ and JF_ o 5_ are continuous bijections between 
X and Xi. The inverse mappings take the form 

(J-±o5±)"' = /+ ^± 

where I is the identity map and 

^±:Xi^L^(M)nC(M) 

are continuous maps. Moreover, both the direct maps T± o S± and their inverses are 
analytic in the sense of Definition \2.1\ below. 
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We note that the mappings of the above theorem are not analytic in the usual sense 
as, for example, the direct maps depend on u and u and the inverse maps depend on F 
and F. 

Now let be the restriction of X to real- valued functions, and let (Xk)-^ be 
the analogous restriction of Xi. If u is real- valued, the reflection coefficients obey the 
additional constraint 

^^ = r±(-s), (1.17) 

which implies that -F+ and -F_ are real- valued. Thus, the maps J^±oS± restrict to maps 
from to (Xjr)-,^. From Theorem 11.11 and the analyticity properties of the direct and 
inverse maps, we obtain our main result on inverse scattering for Schrodinger operators 
with Miura potentials: 

Corollary 1.2 The maps J-'^oS^ and J-'_oS^ restrict to real- analytic bijections between 
Xk and (X^)^. Moreover, 

(J-±o5±)-' = /+ ^± 

where I is the identity map and 

^± : (Xk)i ^Li(M)nC(M). 

Remark 1.3 The class of Miura potentials with a Riccati representative u G X^ is 
"non-generic" among potentials generating non-negative Schrodinger operators: even 
for regular potentials, the generic case is that r(0) = —1 (see, e.g., [13], Theorem 1, 
p. 146 and Remark 9, pp. 152-153, and Example \1.7\ below). In the generic case, the 
scattering problem can be parametrized by left and right Riccati representatives in L^(M) 
having additional integrability on left and right half-lines [24]. 

In the proof of Theorem 11.11 we reconstruct m on a half-line (c, 00) using r_|_ and 
a "right" Gelfand-Levitan-Marchenko equation, and we reconstruct m on a half-line 
(—00, c) using a "left" Gelfand-Levitan-Marchenko equation, and show that these 
reconstructions give continuous maps from Xi into the respective spaces X^ : = 
L^{c, 00) n -L^(c, 00) and X~ := 00, c) fl 00, c). We then exploit the existence 
of an involution X that intertwines 5+ and iS_ to obtain continuity of the inverse maps 
into the full space X. 

Inverse scattering for the AKNS system has been extensively studied: see for 
example the original papers [41] and [1], Shabat's solution of the inverse problem 
for potentials in L^(]R) [38,39], the works of Beals and Coifman [5-7], Cohen and 
Kappeler [11], Wadati [40], and the monographs of Faddeev and Takhtajan [18], and 
Beals, Deift, and Tomei [8]. Aktosun, Klaus, and van der Mee [3] investigated direct and 
inverse scattering for AKNS systems on the line with integrable matrix- valued potentials 
and, in particular, derived partial characterization of the scattering data. Earlier in [2], 
they extended the classical scattering theory for Schrodinger operators by treating a 
class of singular potentials constructed via the Darboux transformations that lead to 
ambiguities in the inverse scattering; see also the related paper by Degasperis and 
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Sabatier [12]. Recently, Demontis and van der Mee [15] studied scattering theory and 
wave operators for AKNS-type equations including those considered here and completely 
characterized scattering data for potentials Q G L^(]R). Novikov [36] studied inverse 
scattering for u in L^-based Sobolev spaces iy"''^(M) with n > 1, and showed that 
the leading singularities of u (modulo iy"~^'^(]R)) can be recovered from the scattering 
coefficients via the Fourier transform. In a similar way, Theorem 11.11 can be interpreted 
as saying that the singularities of u are recovered modulo C(]R) from the Fourier 
transform of the reflection coefficient. Xin Zhou [42] used Riemann-Hilbert techniques 
to show that the direct scattering maps S± : u ^ r± are bijections from H^'^{M) to 
ijf'^(M) where H^^^{R) is the weighted Sobolev space of functions with j derivatives 
belonging to L^(M, (1 + \x\'^)^dx), and 

H{'\R) = {r G H^^\R) : ||r||oo < l} • 

The H^'^ classes are very natural in that certain such classes are preserved under the 
completely integrable flows for the NLS equation, the mKdV equation, etc. On the other 
hand, it is of interest to see how far these can be relaxed and still obtain a bijective 
map. While our class is not preserved under KdV or mKdV, the mappings ^± extend 
to the space 

{L\W))^ = {Fe L\R) : F G L~(M), < l} • 

This space is preserved under maps of the form F i— e''**'*F. In a separate paper [10], 
we will use this extension to construct Strichartz solutions to the defocusing cubic NLS 
equation and study their asymptotic behavior. 

We close this introduction with several examples of Miura potentials. 

Example 1.4 Let u he an even function that for x > equals x~°sinx^. Assume that 
a > 1 and (3 > a + 1. Then u belongs to X and the corresponding Miura potential 
q = u' + is of the form 

q{x) = (3 sign (a;)|a;|^~"~"^ cos \x\'^ + q{x) 

for some hounded function q. Thus q is unbounded and oscillatory; nevertheless, the 
corresponding Schrddinger operator possesses only absolutely continuous spectrum filling 
out the positive semi-axis and the scattering and inverse scattering on such potential is 
well-defined. 

Example 1.5 Assume that G C^(M) is such that (p = 1 on (—1, 1). Take a function 
u{x) = a4>{x)log\x\ with a > 0. Then u G X; moreover, since the distributional 
derivative o/log|a;| is the distribution P. v. 1/x, the corresponding Miura potential q is 
smooth outside the origin and has there a Coulomb-type singularity. See e.g. [9, 19,27] 
and the references therein for discussion and rigorous treatment of Schrddinger operators 
with Coulomb potentials. 

Example 1.6 (Frayer [20]) The Riccati representative u = a;X[-i,i], with a a nonzero 
real constant and xa the indicator function of a set A, corresponds to the Miura potential 

q = a6{- + 1) - a6{- - 1) + 
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6 being the Dirac delta-function centered at the origin. One can explicitly calculate the 
Jost solutions and hence the reflection coefficients by matching solutions in the regions 
X < —1, —1 < X < 1, and x > 1 at x = —1 and x = +1; however, the formula is rather 
involved. We omit details. 

Example 1.7 For the potential q{x) = aS{x), with a a nonzero real constant, the 
extremal solutions ip± of the zero-energy Schrddinger equation (see Section \2.2\) are 
different and equal 

for X > 0, 
for X < 0, 



i^+{x) 




and 



^_(x) 




for X > 0, 
for X < 0, 

respectively (see also [26, Appendix A] and [28]). Therefore such q does not belong to 
the class of singular potentials treated in this paper. 

For this potential, the reflection coefficient satisfies r(0) = —1. This is in fact a 
generic situation even for potentials q in the Faddeev-Marchenko class L^(]R, (l+|x|) dx), 
for which the standard scattering theory is well understood. Therefore a generic Faddeev- 
Marchenko potential is not a Miura potential in the sense of this paper. 

As mentioned above, we shall extend the scattering theory developed here to a wider 
class of potentials including delta functions, Faddeev- Marchenko potentials and many 
other in the papers [24] and [25]. 

This paper is organized as follows. In ^ we introduce necessary notation and 
function spaces and recall some basic results about the Miura map. We review the 
direct scattering problem in ^and prove Theorem II. II in ^ In an appendix, we review 
the Wiener and associated algebras used to analyze the direct and inverse scattering 
maps. 



2. Preliminaries 

2.1. Notation 

In what follows, M2(C) will stand for the 2x2 complex matrices with the Euclidean 
operator norm | ■ | and 




is the usual Pauli matrix. Upright capital Roman or Greek letters (e.g.. A, I, \1/ etc) 
will usually denote elements of M2(C); in particular, I is the 2x2 identity matrix. 

Next, L^(M) and L^(R) are the standard Lebesgue function spaces on the real line M, 
and II ■ 111 and || • ||2 will stand for the respective and L^-norms of scalar, vector, or 
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matrix- valued functions on R. Recalling the Banach space X (see 01.141] and fll.lSp ). 
we define, for any c G M, 

X+ = Li(c, oo) nL2(c,oo) 

and 

X- = L\-oo, c) n L2(-oo, c). 

As usual, C^(]R) denotes the linear space of all functions of compact support that 
are infinitely often differentiable. The Schwartz class S'(M) consists of all infinitely often 
differentiable functions / such that x*^/*-'^ are bounded on M for every natural k and /. 
H^{R) is the space of functions in L^(M) whose distributional derivative also belong 
to L^(M), the norm being 



\m = {\\J ||2 + II/'||2J 

The space H~^{Wj is a dual space to H^{M.), i.e., the space of all distributions that are 
bounded functionals on if^(R). If Y is any of the above spaces, then Y|oc is defined as 
the space of all functions (or distributions) / such that f(j)EY for every G C^(]R). 

Finally, we denote by (M) the Hardy spaces of functions on the real line having 
respective representations of the form 



POO 

Uis)= / e-f<7+(C)dC, 
Jo 

f-{s)= I e-f^?_(C)dC 

J — oo 



for g+ G L2(0,oo) and G L2(-oo,0). Clearly, L'^{M) = Hl{R) ® Hl{R). For 
/ G -^v^(]R), we denote by C the Cauchy integral operator 

(Cf) {z) = ^ [ -^f{s) ds, zeC\R, 
2m J^s- z 

and by C± the operators 

(C±/)(s) = lim(C/)(s±z£), seR, (2.1) 

where the limit is taken in L^(M). The operators C+ and — C_ are orthogonal projections 
onto if^(M) and if^(M) respectively, and C+ — C_ = / as operators on L^(M). We have 
the formulas 

(C±/) (s) = ±^-ix±^+/, (2.2) 

where x+ (resp. X-) is the indicator function of (0, oo) (resp. of (— oo,0)). 

Both the direct and inverse scattering maps are obtained as solutions of integral 
equations that naturally lead to series representations. To capture their analyticity 
properties, we make the following rather nonstandard definition (see [16, 35] and [37, 
Appendix A] for discussion of analytic mappings between Banach spaces). 
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Definition 2.1 Let X he a Banach space of complex-valued functions, let U be an 
open subset of X containing along with every f its complex conjugate f , and let be 
a nonlinear mapping from U into a Banach space Y . We say that ^ is analytic in 
variables f and f if !^(/) = !^2(/, f), where \l/2 : U x U Y is analytic. 

Remark 2.2 If >]/ : X Y is analytic in f and f , then the restriction of 1/ to the 
subspace Xr of real-valued functions is a real-analytic map from X^ to Y . 

We also note that if I' is analytic in / and / and is an analytic function of a 
complex variable, then the point-wise composition (p o fp often turns out to be analytic 
in / and /, see Lemma [A. 3[ 

2.2. The Miura map and Riccati representatives 

In this subsection we briefly recall some results of [26] relating the Miura map and 
the set of positive solutions of the Schrodinger equation at zero energy. We begin 
by considering distributional solutions of the Schrodinger equation with a real-valued 
potential q E H^l(R). For such q, we define the quadratic form s via (11.21) and let 



The quadratic form s is nonnegative if Xo{q) > 0. A function y G is a 

{distributional) solution of —y" + = if for every ip G C^(]R), 3{y,ip) = 0, with 
s{y, if) denoting the value of the associated sesquilinear form. It is not difficult to see 
that any solution y may only have isolated zeros where the solution y changes sign; in 
particular, if y is nonnegative, it is either identically zero or strictly positive. If Ao(g) > 
and ?/ is a solution with either y G L^(M+) or ?/ G L^(M~), then y has a single sign or 
is identically zero: if y has a zero, then there is a half-line Dirichlet problem operator 
—d?/dx^ + q with a negative eigenvalue, contradicting the assumption that Ao(g) > 0. 
Finally, \o{q) > if and only if the set 



of positive distributional solutions to the equation —i})" + g?/; = is nonempty. If 
ip G Pos (g), then q = u' + v? ioi u = ip' 

We initially consider the Miura map B : u ^ u' + v? acting from Lj^q^(]R) 
into ifj^^(M), by interpreting u' as the distribution derivative of u. If q lies in the 
range Rani? of -B, then any function u G B^^{q) is called a Riccati representative of q. 
The range of the Miura map consists of those q G iJ[~^(M), for which Ao(g) > 0, or, 
equivalently, for which the set Pos (g) is nonempty. The set Ran B is nowhere dense 



Ao(g) :=inf{s(¥p) : ^eC^iR), M2 = l}. 



Pos(g) = {ije Hl^M) : ^ > 0, m = 1, Vv^ G C^iR) 5{^,ij) = O] 





and 




Inverse scattering for Miura potentials 



11 



take B~^{q) into Pos (g) and vice versa, so that these two sets are homeomorphic 
respectively as subsets of Lj^qj,(M) and Hl^^{M.). Given any yi G Pos(g), any solution 
of —y" + qy = with y{0) = 1 (whether positive or not) takes the form 

y{x) = yi{x) + c yi{sy'^ ds^ 

for some c (see [22, Ch. IX. 2 (ix)]). Using this fact, it is not difficult to see that the set 
Pos (g) is the set of convex combinations of extremal solutions and (f- characterized 
by the respective divergence of f^(p^'^{s)ds and <^I^(s) ds. As we will see, the 
extremal solutions (p+ and ip- coincide up to constant factors with the usual right and 
left Jost solutions at zero energy whenever the latter are defined. The extremal solutions 
ip^ and if - either coincide or are linearly independent, so that B~^{q) (with the topology 
induced from Lj^Qp(M)) is homeomorphic either to a point or to a line segment. 

Here we are concerned with the restriction of B to L^(]R) so that the range of B is 
contained in H~^{W). A distribution q G H~^{M.) lies in the range of B if and only if 
the positivity condition Xo{q) > holds and q can be presented as f + g for / G L^(M) 
and g G L^(M). The following result is a consequence of the proof of Lemma 4.1 in [26]. 

Lemma 2.3 Suppose that q G H^^{Wj and that there exists a Riccati representative u 
of q with u G L2(]R+). Then the same is true for every Riccati representative of q. A 
similar statement holds with //^(M"*") replaced by L^(R^). 

The following simple uniqueness result is very important for the present paper. 

Lemma 2.4 Suppose that q G H^^(R) and that ui and U2 are Riccati representatives 
of q belonging to Xq . Then Ui = U2- Similarly, if Ui and U2 are Riccati representatives 
of q belonging to Xq , then Ui = U2- 

Proof If u G X^ = L\R+) f] L2(M+) belongs to B'^q), then (p{x) := exp {J^ u{s) ds) 
is bounded above and below by fixed constants for a; > 0, so ip{s)~'^ ds diverges. 
Since ip solves the equation —y" + qy = 0, we conclude that ip = ip+ is an extremal 
positive solution, and u = Lp'^jip^^. Hence u\ = U2 = ip\/ip+. The other proof is similar. 

□ 

Recall that X = L^(R) fl L^(]R). We immediately obtain from Lemma [23] the 
following uniqueness result: 

Corollary 2.5 Every q G H^^{M.) has at most one Riccati representative u belonging 
to X. 

3. Direct Scattering 

Although direct scattering for the ZS-AKNS system is well-understood (see, for example, 
Beals and Coifman [5], Zhou [42], and the monographs of Faddeev and Takhtajan [18] 
and Beals, Deift, and Tomei [8]), we give a brief synopsis here. We will study the direct 
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scattering map, derive representation formulas for the scattering solutions, and finally 
derive the Gelfand-Levitan-Marchenko equations and the reconstruction formula. 

It is useful to factor out the leading behavior of solutions to (11. 3p with Q of (11.111) 
by setting 

\['(x, z) = M(x, z) exp(iza;(j), 
where M takes values in 2 x 2 complex matrices, so that 

M; = izad,(M) + Q(x)M, (3.1) 

where 

adA(B) := [A, B] = AB - BA. 
To study scattering asymptotics, we remove the linear term from (13.11) by setting 

M(x, z) = exp(i2;x ado-) N(x, z), (3.2) 

where 



exp(it ado-) 



a h \ I a e**6 
c d J \ e~**c d 



is an automorphism of the Lie group 

SU{1, 1) := {Y e M2(C) : YVY = a, det Y = 1} . 
Then N(x, z) obeys the differential equation 

^'^{x,z) = Q{x,z)^{x,z) (3.3) 

where 

(0 e'^'^'^uix) \ 
Q u\x) U / 

Since G is traceless and crG + Gcr = 0, it is not difficult to see that det N and N*crN are 
independent of x for any solution of (13.31) . 

Now consider the singular initial value problem 

N^(x,z)=G(x,^)N(x,z) 
N(oo,2) =1. 

Reformulating this problem as an integral equation 

N(x,z) = l+ / G(y,z)N(|/,z)dy, (3.6) 



it is not difficult to see that if m G (M) and z = sv& real, there exists a unique solution 
which is a continuous curve in S\J(\^ 1) with a well-defined limit as x — — oo. Tracing 
through the definitions and comparing with (11.61) . we see that 

A(s) = N(-oo,s). 
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Since N(— oo,s) G SU{1, 1), we may write 

N(-co..)=f (3.7) 

In what follows, we will study the direct scattering map by studying the singular initial 
value problem (13.51) . We will then return to (13. ip in order to study the scattering 
solutions and their analyticity properties. From these properties we can deduce the 
Gelfand-Levitan-Marchenko equations and reconstruction formulas for u. 



3.1. The Direct Scattering Map 

The basic objects of direct scattering are the reflection coefficients (I1.7P and (11.81) 
computed from the asymptotic behavior of solutions to (11.31) . It is simplest to study the 
map 



u 



a{s) b{s) 
b{s) a{s) 



using (13.51) and (13.71) and deduce properties of the reflection coefficients from those of a 
and b. 

The first result is: 

Lemma 3.1 Suppose that -u G X. Then the representations 

POO 

a(s) = 1+ / e'^^A(C)dC, (3.8) 
Jo 

/oo 
e'<B{C) dC (3.9) 
-oo 

hold, where A,B E X with 

Pill, ||5||i<exp(||^i||i)-l, (3.10) 

P||2,||5||2 <exp(||^i||i)||^i||2. (3.11) 

Moreover, the maps u ^ A and u ^ B are analytic in u and u in the sense of Defini- 
tion\2l\ 

Proof. To obtain the representations (I3.8p - (l3.9p and the estimates (I3.10p - (l3.1ip . we 
study the Volterra series for (13. 6p . A straightforward computation with (13.60 using 
(13. 4p shows that the entries rtn and n2i of N obey the system of integral equations 



nii(x, z) = l+ / e-''yu{y)n2i{y, z) dy, (3.12) 

J oo 

n2i{x,z)= / e'^%(?/)?2ii(?/, z) dy. (3.13) 



We can iterate equations (I3.12p - (l3.13p to obtain Volterra series representations, and 
then obtain the required formulas for a and b since, by (13. 7p . we have a = nii(— oo, • ) 
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and b — 7121 (—00, • ). First, we have 

00 

nu{x,s) = l + '^T2n{x,s), (3.14) 

n=l 

00 

n2i{x,s) = ^T2n+i{x,s). (3.15) 

n=0 

Here 

Tn{x,s):^ ex.p{izipn{y))Un{y) dy, 

J x<yi<---<y n 

where y := (|/i, ...,yn), 

n-l 
3=0 

U2n{y) ■= U{yi)u{y2) ■ ■ ■U{y2n-l)u{y2n): 

and 

U2n+i{y) — -u{yi)u{y2)u{y2) • ■ ■ u{y2n)u{y2n+i) ■ 

Note that <^2n(y) > if |/i < • • • < y2n, whereas (/?2n-i(y) > x for y in the range of 
integration for T2n-i- We may write 

POO 

T2n{x,s) = / e'^«^2n(2:,C)dC, 

T2n+l{X-, s) 

where 

A2„(X,C) := / C/2n(y)d52n, 

a:<S/l<---<y2n 

B2n+l{x,C)-= U2n+l{y) dS2n+l 

J 



POO 

/ e'<B2n+iix,C)dC, 

J X 



V2n+l(y)=f 
x<j/i<---<a2n+l 



are multihnear functions of u and u. Here djS'n is the surface measure on the 
hypersurface {y e IR" : </?„(y) = C}- Using the fact that dy = d-S^dC and letting 77 
and 7 be the bounded monotone functions 



i^{x) := / \u{t)\dt, 

1/2 



X 



7(x) := (^£\u{t)\Ml^ 
we can easily derive the estimates 

l|B,.«(.,.)lk<|^. 
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\\A2n{x.-)h <^^7(^), 
l|i?2„+l(x, •)I|2<^7(X). 

These estimates show that the Volterra series f l3.14l) - p.l5l) converge and define 
bounded absolutely continuous functions of x having finite limits as x — — oo for each s. 
By taking hmits in fl3.14p - fl3.15p . it is then easy to deduce that a and h have the series 
representations 

oo 

a{s) = 1 + y^T2n(g), 

n=l 

oo 

^(*) = X^^2n-l(s), 



where 



n=l 



T„(s) := lim T„(a;, s) 

X— oo 



is a bounded continuous function of s for each n. Setting A2n{C) = A2n{-~oo,Q and 
52„_i(C) = 52„_i(-oo,C) and defining 

oo 

n=l 

oo 

5(C):=5^i?2n+l(C), 

n=0 

we obtain the required representations. The continuity follows from the estimates 

||^2n-^2n||p < 1- TT^Wu - u\\p, 

[2n — 1)! 

~ C^" 

||-B2„+i — -B2n+i||p < tttttII^ ~ ^llp 
(in)! 

for p = 1,2 where C := HmHi + □ 

Since A(s) G ^[/(l, 1) for every s G M, we have the identity 

\a{s)\'-\b{s)\' = l, (3.16) 

and thus a does not vanish on the real line. Moreover, the integral representation (13.80 
shows that the function a is the boundary value of a function of the complex variable z 
(also denoted by a) that is defined in the closed upper half-plane C+ by the formula 



POO 

a{z) :=1+ / e^<A(C)dC. 
Jo 



The function a so defined is continuous in C+, analytic in C"*", and has no zeros there 
(see, e.g., [1], [14, Sect. 3]). 
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Remark 3.2 Since a does not vanish on the real line and tends to 1 at infinity in view 
of ( 13*. gj) and the Riemann-Lebesgue lemma, it follows from Lemma \A.^ that a is an 



invertible element of the Banach algebra l+X (see Appendix A ). The function can 
be extended to the upper half-plane as an analytic function l/a{z); thus belongs 
to the Hardy space Hi and 



POO 

\z) = l+ / C(C)e'<dC 

^0 



for a function C G . Moreover, C is analytic in u and u in the sense of Definition \2.1\ 
as follows from Lemma \A.3\ by setting (j){z) := 1/z. 

Recalling the formulas (11.71) and (ll.Sp . the constraint |a(s)|^ — = 1, and the 

bound ||a||oo < 1 + ll^lli) it is easy to see that the reflection coefficients r_|_ and r_ 
are continuous with ||r-|-||oo < 1- By Remark I3.2[ r± are also elements of the Banach 
algebra X (i.e., Fourier transforms of functions in X) which are analytic in u and u. 
Thus: 

Lemma 3.3 The maps J^^oS^ and J^^oS^ are continuous maps of X into Xi, analytic 
in the sense of Definition \2.1[ 

It will be important to formulate the exact relationship between r+ and r_. It 
follows from the definitions (ll.yp -l ITISl) that 



-r+(sj 



a[s] 



a{s) 

The relation (13.161) yields for s G M the identity 

\a{s)r = 1 - \r{s)\^ 

where r = r_ or r = r+, so that \a\ is determined by r± on the real line M. Since the 
function a is bounded and continuous in C+, analytic in C"*", and has no zeros in C+, 
the function log a is well defined in C+ and analytic in C^. Observe that 

Re loga(s) = log \a{s)\ = ~ log(l - \r{s)\^); 

thus we can use the Schwarz formula to reconstruct the function log a from its real part 
on M. Explicitly, we get 

r 1 /"°° 1 1 
a{z) = exp — / log (l - |r(C)|^) dC , 

where r = r+ or r_; the boundary value is 

a(.)=exp{C+[log(l-|r(-)P)](5)}, (3.17) 

where C+ is the Cauchy operator (12. ip . which is bounded from L^(M) to itself. We also 
observe that the problem of reconstruction of the function a in C"*" from the values of \a\ 
on the real line could be reduced to the classical Riemann-Hilbert problem, cf. [14,42]. 

Let us denote by Xi the image of Xi under the Fourier transform (I1.16p . If r G Xi, 
it follows from Banach algebra properties of X (see Appendix A ) and Lemma [A. 21 that. 
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also, log (l — |?"(C)|^) £ X. By Lemma [A.31 log (l — |?"(C)|^) is analytic in r and r in 
the sense of Definition I2.1[ From fl2.2p we deduce that the function 

g{s)=C^ [log(l-IK-)r)] {s) 

also belongs to X and is analytic in r and f in the sense of Definition 12. 1[ It now follows 
from Lemma [A. 2 1 that 

a{s) - 1 = exp[5((s)] - 1 

belongs to X as well and is also analytic in the same sense. Hence, a of fl3.17p is an 
invertible element of 1 + X whenever r E Xi. 
Now we define a nonlinear map 

I:Xi^Xi 

by 

(Xr)(.) = -K^ 

a[s) 

where a is given by (13 . 1 7p . Note that in the Schrodinger case, where (11.170 holds and, 
also, a{—s) = a{s), the map X is given by 

, a(—s) 



a{s) 

(cf. for example [13], §2.2, Theorem 1). 

Lemma 3.4 The map I is a continuous involution of Xi, analytic in r and f in the 
sense of Definition \2. l[ 

Proof. The map r y-^ a defined by the solution of the Schwarz problem above is 
continuous from Xi to 1 + X and a ^— a/a is continuous from 1 -j- X to itself. 
Multiplication by r is continuous from 1 -j-X to X, so we conclude that X is continuous 
as claimed. To see that X is an involution, we note that |(Xr)(s)| = |r(s)| and deduce 
from ( ]3.17p that X^r = r. Finally, the analyticity follows from analyticity of the map 
r ^ a and invertibility of a in 1 -j- X. □ 

It now follows that the identities 

S-=IoS+ 

and 

5+ = X o 5_ 

hold. 
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3.2. Scattering Solutions 

Recalling the scattering solutions \E'± of ( 11. 3p with 

lim \^±{x,z) — exp(i2;xcr)| = 0, 

X— >±oo 

let us write 

\E'±(x, z) = M-|-(a;, z) exp{izxa), 
so that M± solve (13.11) with 

lim \M±{x,z) - I| = 0. 

X— utoo 

These singular initial value problems can be reformulated as Volterra integral equations 
M±(x, z)=l+ I exp(i2;(a; - y) ad,) [Q){y)M{y, z)] dy. (3.18) 

J±co 

These equations have a unique absolutely continuous solution for each real z provided 
that u G L^(M). The map x i— >■ y\.±{x,z) is a continuous curve in SU{1,1) for each 
fixed z. Letting mf- denote the entries of M-t, it is not difficult to see that 

mf^{x,z) = l+ u{y)m^^{y,z)dy, (3.19) 



m+(x,^)= / e-'<''-yMy)mUy,z)dy (3.20) 
(and similar formulas with all + signs replaced by — signs), while 



"^12(3^5 z) = mj^{x, z), (3.21) 

^22 (2^) = ^iiix, z) (3.22) 

(and similar formulas with + replaced by — ) since M-t take values in SU{1, 1). Thus 
the matrix-valued function M_|. is determined by solutions of the system of two integral 
equations f l3.19l) - p.20l) . and M_ is determined by solutions of a similar system for m'^^ 
and m^]^. 

The integral equations fl3.19p - fl3.20p can be analyzed in the same way as equations 
f i3.i2p - fi3.i3p in the proof of Lemma [3. 11 or we can note from (13. 2p that 

m'l^{x^z)= nj,^(x,2;), (3.23) 

m+ (x, z) = e-'^^n+ (x, z). (3.24) 

It follows directly from the representations fl3.14p - fl3.15p already established for nfi and 



that 



/■oo 

m+(x,^) = l+ / e'^^r+(x,C)dC, 

poo 

mt,{x,z)= / e'<r+(x,C)dC, 
^0 



where the estimates 

||rn(x, Olli <exp [ri{x)] - 1, 
||rn(x, ■ )ll2 < exp ['r]{x)]-f{x) 
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hold (note that in the series for mjj^, the factor exp(— i2;x) in fl3.24p changes the phase 
function to v'2n+i(y) — this phase function is nonnegative for any x). 

The exphcit multihnear expressions for T^-^ (x, ■ ) and {x, ■ ) show that they 
are analytic functions of u and u in the sense of Definition 12.11 taking values in the 
space Xq for each real x. Using the symmetry properties f l3.2ip - p.22p . we recover 
the full matrix- valued function r+. We can make an analogous analysis for M_, and 
obtain the following representations. In what follows, (8> M2(C) denotes the space of 
matrix- valued functions over with entries belonging to X^. 

Proposition 3.5 Suppose that u E X and that M± are the unique solutions of the 
integral equations \3.18\) . Then there exist continuous functions 



^3x^V 



{x,- )e X+ ® M2(C) 



and 



M 3 X ^ r_(x, ■)ex^ ® M2(C) 



with 




and 




so that the representations 




(3.25) 



(3.26) 



hold. Moreover, the estimates 




)||i < exp[r7+(x)] - 1, 
)||i < exp[r7_(a;)] - 1, 
)||2 < exp[r/+(x)]7+(x), 
)||2 < exp[r7_(a;)]7_(a;) 



hold, where 




±00 



Finally, for every fixed a; G M, the maps u ^ T±{x, ■) from X into X^ ® M2(C) are 



analytic in the sense of Definition \2.1\ 



Inverse scattering for Miura potentials 



20 



Remark 3.6 The function r+(x, C) depends only on values of u{y) with y > x, while 
the function r_(x, Q depends only on values of u{y) with y < x. We will derive integral 
equations for T± given the functions F± of l\l.l!!^) - ^1.13\) : the equation for will provide 
a stable reconstruction of u{x) on any half-line (c, oo), while the equation for r_ will 
provide a stable reconstruction for any half-line (— oo,c). 

If we write M_|_ = (mf,m2), where mf and are column vectors, and similarly 
write M_ = (m]^,m^), then (mf,m2) extends to an analytic matrix-valued function 
of z in Imz > 0, while (m^,m^) extends to an analytic matrix-valued function of z in 
Im2; < 0. Indeed, as an immediate corollary of Proposition 13.51 we have: 

Proposition 3.7 For every fixed x G M, the functions 



belong to the Hardy space H^{R,C'^). 

These Hardy space properties are the starting point for the derivation of the 
Gelfand-Levitan-Marchenko equations. 

3.3. The Gelfand-Levitan-Marchenko equations 

We now derive integral equations for F-i- in terms of F± and also derive a reconstruction 
formula for u from r±. In the next section we will solve these equations and prove 
consistency of the two reconstructions. Given functions F± of fll.l2l) - fll.l3l) . let us 
define 



Proposition 3.8 Suppose that u G X; then the kernels and r_ obey the equations 



T_{x,C) + ^-{x + C)+ T_{x,t)n_{x + t + C)dt = for a.e. C<0, (3.27) 



T+{x,C) + ^+{x + + T+{x,t)n+{x + t + C)dt =0 for a.e. C > 0. (3.28) 



The equations ( I3.27P and (13.28^ are the left and right Gelfand-Levitan-Marchenko 
equations. 

Proof. We will derive f l3.27p : the derivation of fl3.28p is similar and will be omitted. We 
begin with the identity 





We will prove: 
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for a fixed x and real z, which implies that 



mi 



r_e '^^m 



2 ) 



mo + r_ e'^^m. 



(3.29) 
(3.30) 



From Proposition 13.71 and Remark 13.21 the left-hand side of (13.291) satisfies 

mt 



1 




for each x. On the other hand, it follows from Proposition 13.51 and ( 11.13^ that 



+ r_e ^^^m2 



1 




is represented as 



r-(x,c) + F-(x + C) 





1 



e^<dC 



/ T~{x,t)F_{x + t + Odt 

'OO J ~oo 



(3.31) 



e^^^ dC, 



where we regard the columns T^{x, ■) and (x, ■) of F_( vector-valued 
functions on the line which vanish for > 0. Since (13.311) should give a function 
in Hi 



^2 ^R,C2), we get 



Fr(x,C) + F_(a; + C) 




1 



F2(x,t)F_(x + t + C)dt = (3.32) 



for almost every < 0. A similar calculation with (I3.30p shows that 

F2-(x,C) + F-(x + C) (| J ) + / rr(x,t)F_(x + t + C) dt = (3.33) 

for almost every ^ < 0. The equations fl3.32p and f l3.33p together imply fl3.27p . □ 

Finally we recall the formulas for reconstructing u from which play a crucial 
role in the inverse theory. 

Proposition 3.9 Suppose that u G ^(M). Then 

m(x) = -F+(a;,0) (3.34) 

and 

m(x) = F72(x,0). (3.35) 

Proof. If u G 5(1^), it is easy to see that the kernels T±{x,C,) are smooth functions of 
X and and that the derivatives are integrable in C,. Using the integral representations 
for M-t in terms of T± together with the integration by parts formula 

k 

/(C) 



exp(izCa3)dC = [-^^^) f'-'K^) + O (^) , 
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true for / G C'=+i([0, oo); M2(C)) with ^ \f^^\C)\^C < oo for < j < + it is easy to 
see from the integral representations (I3.25l) - fl3.26l) that M+ and M_ have differentiable 
asymptotic expansions of the form 

M,(x,^)~I+M±iM + o(l) (3.36) 

for u G 5'(]R), where 

M±,i(x) = ±ir±(x,0)or3. 

Substituting (13.361) into (13. ip . we conclude that 
Q(x) = -iad^(M±,i(x)) 

= ±i[a3r±(x,0)a3-r±(x,0)]. 
The formulas (13.341) and (13.351) are an immediate consequence. □ 



4. Inverse Scattering 

We now suppose given (F+,F_) G Xi x Xi with F_ = JF_(TF+) and show that 
the equations (I3.27P and (13.280 admit unique solutions V± depending continuously on 
F± G Xi. We then argue by continuity and density that the reconstruction formulas 
(I3.34P and (13.350 . known from the classical theory for F± in the Schwartz class 5'(M), 
continue to hold also for F± G Xi. We remark that equations (I3.27P and (I3.28P can 
also be analyzed using the factorization theory in operator algebras as developed e.g. 
in [21, Ch. IV]; see also [33,34] for some extensions and [4] for an example of application 
in inverse problems for Dirac operators in the AKNS form. 

We give the details for recovering u from since the recovery of u from F_ is 
closely analogous. By (I3.34p . it will suffice to solve (13.281) for ^^^d show that the 
solution is sufficiently regular so that u{x) = —Tf2{x,0) is well-defined and gives a 
continuous map from Xi to X+ for any c G M. 

First we prove the existence of a unique solution. From (I3.28p we have (writing F 
for F+) 

POO 

F{x + c) + r+ (x, c) + / r+ (x, t)F{x + 1 + c) dt = 0, 

Jo 

POO 

r+(x,c)+ / r+(x,t)F(x + t + c)dt = o. 

Jo 

Iterating these equations and writing F for F^2) we see that 

POO POO 

T{x,C)+F{x + C)- / T{x,t2)F{x + t2 + ti)F{x + ti + C)dt2dti = 0. (4.1) 

Jo Jo 

To solve this equation, we consider the scalar operators Tp{x) and T-p{x) on L^(M+) 
given by 



[T^(x)V^](C)= / ^(t)F(x + t + C)dt, 
Jo 

POO 
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and regard fl4.1l) as an equation in L^(R"'") for each fixed x: 

r(x, ■) + F{x+ ■)- {Tf{x) o Tp{x)) r(x, ■ ) = 0. (4.2) 
The operators Tf{x) and T-p{x) satisfy the estimates 

\\Tf{x)\\li^L1, \\T-p{x)\\li^L^ < ||F||i, 

\\Tf{x)\\l^^l^, \\T-p{x)\\l^^l'^ < (4.3) 

\\Tf{x)\\l2^L2, \\Tp{x)\\L2^L^ < ||r||oo < 1, 

where r = F. Moreover, the map 
F^Tf 

is a continuous map from X into the bounded operators from L^(M+) to itself, from 
L^{R+) to L2(]R+), and from L^{R+) to itself. Since Tp is compact as an operator from 
L^(]R"'') to itself for F G C^(]R) and p = 1, 2, it follows by density that the same is true 
for any F G X. 

The estimates imply that the operator 

{I-TF{x)oTp{x)y' 

exists as a bounded operator from L^(R"^) to itself given by a convergent Neumann 
series. In particular, keYF2(M.+) {I — Tf{x) oT-p(x)) is trivial. On the other hand, any 
solution tp G L"'^(]R+) of the equation 

{I -TF{x)oTp{x))ij = Q 

actually belongs to L'^{R+) since Tf{x)oTp{x) maps L^{R^) into L2(R+). It follows that 
kerii(R+) (/ — Tf{x) o T-p(x)) is also trivial. Since Tp and Tp are compact, we conclude 
that (J — Tf{x) o T-p[x))~^ also exists as a bounded operator from L-^(M+) to L^(]R+). 
The solution to (14.21) is given by 

r(x, ■) = -{!- Tf{x) o Tp{x))-' F, =: -F, - G{x, ■ ), (4.4) 

where Fx{ ■ ) := F{x + ■ ) and 

G{x, ■ ) := Tf{x) (/ - Tp{x) o Tpix))-' [Tp{x)F,] . 

The right-hand side of (14. 4p defines for each x G M a function r(x, ■ ) G X^. 
To study continuity in x and (, set 

H{x, ■):={I- Tp{x) o TF{x)r' [T^(x)Fj . (4.5) 

The estimate 

\\Hix,-)\\L^M+)<il-\\r\\lr'\\r\U\F\\2 

holds, and x ^ H{x, ■ ) is a continuous mapping from R into L^(R'^), as follows from 
the continuity of Tf{x) as an operator- valued function of x, the continuity of x ^ F^ 
as a mapping from R into L^(R+), and the uniform bounds on the resolvent. From the 
formula 



POO 

G{x,C)= H{x,t)F{x + t + C)dt, 
Jo 
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the Schwarz inequality, and the continuity of H{x, ■ ) as a mapping into L^(M"'"), we 
deduce that G{x^ Q is jointly continuous in x and C and is uniformly bounded. Thus: 

Lemma 4.1 Suppose that F G Xi. For each x G M, there exists a unique solu- 
tion r(x, ■ ) of belonging to Xq . Moreover, 

V{xX) = -F{x + ()-G{x,0, 

where G is a hounded, jointly continuous function of x and (. 

We can now obtain a limit of T{x, () a.s( and compute a putative reconstruction 
:= -r{x, 0) = F{x) + w{x), (4.6) 

where 

/•oo 

w{x) = / H{x,t)F{x + t)dt. (4.7) 



Recall that = L^{c, oo) nL^(c, oo). We will now show that, for any fixed c G I 
the map 

g+ : Xi3F^u+e X+ 
defined by (14.61) is bounded and continuous. We will also show that 

oo 

SAF)=F + Y,SL^iiF,F) (4. 

n=l 

for multilinear maps : X" — > X^. Indeed, set 

H^nix, ■ ) := {T-p{x) o TFix)^-' [T-^{x)F,] 

and 



/•CO 

{F,F)= / F{x + t)H2n{x,t)dt. 
Jo 



Then the claimed representation follows from fl4.5p . (14.61) . and (14. 7p provided that the 
series (14.81) converges in 

To prove convergence, suppose that F G Xi and choose p < 1 such that ||-F||oo < P- 
First, from (14.31) we have, 

||i^2n(x, ■)l|2<p'"-ii^||2, 

SO that by the Schwarz inequality 

\\GUiiF,F)\\oo < p'^-'llFWl (4.9) 

It follows that G2n+iiF, F) converges in L°°{R). On the other hand, we have the 

explicit formula 

OL+iiF,F){x) = [ F{x + t,)F{x + ti + ta) ■ • •F(x + t2„-l + t2n)F{x + t^n) dt, 
with t := (ti, . . . , t2n) G Fix F and choose xq so that 

/>oo 

/ \F{s)\ds<p. 

Jxo 
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Note that the same condition holds for some relatively open set of F G Xi. For any F 
in such an open set we then have 

II^2ViIUh-o.oo)<p'"+^ (4.10) 

Thus, combining (14.91) and fl4.10p . we have for any c G M the estimate 

\\Qtn+i\W{c.oo) < (max{xo-c,0}||F||2 + pV""\ (4.11) 

which shows that Xl^i ^2n+i(-^5 converges in L^{c, oo) for any fixed c G M. It now 
follows from (14. 9p and (14.111) that the series (14.81) converges in L^(c, oo), and hence in X+ 
as claimed. 

We have shown: 

Lemma 4.2 For any fixed c G M, the mapping Qj^ : Xi is continuous and admits 

the representation d^.gp . Moreover, — I is a continuous map of Xi into C{M). 

Similarly, we can study the "left" Gelfand-Levitan-Marchenko equation (13.271) to 
define a map Q_ : F_ via the reconstruction formula (I3.35p . The analogue of (14.11) 

(where now F denotes ^"^^ P denotes F_) is 

/ F(x, t2)F{x + t2 + ti)F{x + ti + C) dt2 dti = 0, 

-oo J —oo 

and we set 

M_(x) = limF(x, y) 

where the limit is taken in X~ . In an analogous way we obtain a representation 

oo 

= g^{F) := -F + Y^g^n+AF^F). (4.12) 

n=l 

where Q~ are multilinear functions acting from (Xi)" to X~ . An analogous argument 
shows: 

Lemma 4.3 For any fixed c G M, the mapping ^_ : Xi — > X~ is continuous and admits 
the representation lj^4-12\ ). Moreover, — I is a continuous map of Xi into C(M). 

With these continuity statements, we can prove Theorem II. 1[ 

Proof of Theorem 11.11 We will discuss only the mapping JF_,_ o Sj^. To show it is 
injective, assume that -F+ G Xi corresponds to m G X. Then, solving the Marchenko 
equation (13.280 . we find the kernel F+ that by (I3.25P determines uniquely the Jost 
solution which, in turn, gives the matrix potential Q of (II. lip . Therefore different 
u are mapped into different F+, so that there is at most one potential m G X with a 
given right reflection coefficient. 

We next construct the inverse of o 5+. Given F in Xi, we set u+ = Q+F and 
U- = Q-{J-'-{TF)). We need to show the following facts: 

(1) = and 

(2) the potential « = = has reflection coefficients r+ = F and r_ = XF. 
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Then continuity of the mapping F \^ u follows from those of Q± and (1). Analyticity 
of this map follows similarly from (1), analyticity of ^_ on and ^+ on Xl^^, and the 
representation u = (u^ + (1 — Quj^ for a function ( G with ((t) = 1 for t < — 1 
and C(^) = for if: > 1. Fact (2) shows that the range of J-'+ o 5+ is Xi. 

We will appeal to standard results for the Gelfand-Levitan-Marchenko equations 
when F G S(R) (see for example Section II. 4 of [18]) and the continuity of the direct 
and inverse scattering maps established above. Suppose given F G Xi and let be a 
sequence from S'(M) with F„ — F in X. Note that, for n sufficiently large, F„ G Xi since 
F„ ^ F in Let F+ := F, define F_ as T-{TF), and similarly let F„,+ := F„ 

and define F^ _ as JF_(XF„). 

Now set := ^±F„-|-. By standard theory, we have u±^n ^ 5(1^) with 

= M-,n = Mn, and the potential m„ has associated refiection coefficients Fn^+[s) 
and F„^_(— s). By Lemmas 14.21 and 14.31 Un,+ m+ in X+ and «„_ m_ in X~, so 
that = = u and m„ — > -u in X. By the continuity of the direct scattering maps 
of Lemma 13731 we can also conclude that u has refiection coefficients r+(s) = F+(s) and 
r_(s) = F„(-s) = (JF)(s). □ 

Now it is easy to deduce Corollary 11.21 

Proof of Corollary 11.21 We only need to prove that the maps are real-analytic. Real 
analyticity of the direct maps follows immediately from Lemma 13.31 and Remark 12.21 
Real analyticity of the inverse map follows similarly. □ 
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Appendix A. 

Let us denote by L^(R) the Wiener algebra of Fourier transforms (11.161) of functions 
in L^(M) with norm ||/||;^ := and by X the Banach algebra that is the 

image of X = L^(]R) fl L^(M) under the Fourier transform, equipped with the norm 
WfWx ■~ 11/11^- ^l^o denote by 1 -fX the unital extension of X obtained by adding 
the constant functions and norming 1 -j- X with the norm 

||c + /IL:^ = |c| + ||/||^; 
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we similarly define 1 -j- L^(R). The Fourier transform extends to 1 -j- X by mapping the 
constant 1 into the convolution identity 5. 
We will need the following results. 

Lemma A.l Suppose that / = a + /;,Gl-i-X. Then f is invertible in the Banach 
algebra 1 + X if and only if f is non-vanishing on M and a ^ 0. 

Proof. If / is invertible in 1 -(- X it is also invertible in 1 -j- L^(M), so the condition is 
necessary by the Wiener theorem. If, on the other hand, / does not vanish on R and 
a 7^ 0, then / is invertible in 1 -f L^(M) with /"^ = a"^ + g for g E L^{M.). We need to 
check that g G L^(M). Without loss we take a = 1 and compute that ^ = —(1 + h)~^h, 
which shows that g G L^(R) as required. □ 

We now have an analogue of the Wiener-Levi theorem for 1 -\- X. 

Lemma A. 2 Assume that / G 1 + X and that (j) is a function that is analytic in an 
open neighborhood Q of the closure of the range of f . Then (pof g l + X and, moreover, 
the map 

f ^<P° f 

is an analytic map from 1+X into itself when restricted to functions with range contained 
in a fixed compact subset of Q. 

Proof. It suffices to note that, according to the above, the spectrum of / in 1 -j- X 
coincides with the closure of its range. Then the standard functional calculus for Banach 
algebras applies, thus yielding the result. □ 

Lemma A. 3 Assume that U is an open subset of 1 + X and that I^ : U —>■ 1 + X 
is analytic in f and f in the sense of Definition 2.1. Let also (j) be a complex-valued 
function that is analytic in a neighbourhood of the closure of the set IJ^^^ Ran (/^(Z). 
Then the map / t— > o I/{f) from U into 1 \- X is analytic in f and f . 

Proof. By Lemma IA.2t t-^ o (7 is an analytic mapping in the Banach algebra 1 + X 
defined on I^{U). It is therefore analytic in 1 -j- X considered as a Banach space. Thus 
the mapping f v-^ (j)o I^{f) is analytic in / and / as a composition of analytic mappings 
between Banach spaces. □ 
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